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ABSTRACT 

Geodetic control is extremely important in the production 
and quality control of topographic data sets, enabling 
elevation results to be referenced to an absolute vertical 
datum. Global topographic data with improved geodetic 
accuracy achieved using global Ground Control Point 
(GCP) databases enable more accurate characterization of 
land topography and its change related to solid Earth 
processes, natural hazards and climate change. The 
multiple-beam lidar instrument that will be part of the 
NASA Deformation, Ecosystem Structure and Dynamics of 
Ice (DESDynl) mission will provide a comprehensive, 
global data set that can be used for geodetic control 
purposes. Here we illustrate that potential using data 
acquired by NASA’s Ice, Cloud and land Elevation Satellite 
(ICEsat) that has acquired single-beam, globally distributed 
laser altimeter profiles (± 86°) since February of 2003 [1,2]. 
The profiles provide a consistently referenced elevation data 
set with unprecedented accuracy and quantified 
measurement errors that can be used to generate GCPs with 
sub-decimeter vertical accuracy and better than 10 m 
horizontal accuracy. Like the planned capability for 
DESDynl, ICESat records a waveform that is the elevation 
distribution of energy reflected within the laser footprint 
from vegetation, where present, and the ground where 
illuminated through gaps in any vegetation cover [3]. The 
waveform enables assessment of Digital Elevation Models 
(DEMs) with respect to the highest, centroid, and lowest 
elevations observed by ICESat and in some cases with 
respect to the ground identified beneath vegetation cover. 
Using the ICESat altimetry data we are developing a 
comprehensive database of consistent, global, geodetic 
ground control that will enhance the quality of a variety of 
regional to global DEMs. Here we illustrate the accuracy 
assessment of the Shuttle Radar Topography Mission 
(SRTM) DEM produced for Australia, documenting 
spatially varying elevation biases of several meters in 
magnitude. 

Index Terms — ICESat, global geodetic control, laser 
altimetry, DEM accuracy, elevation errors, DESDynl 


1. INTRODUCTION 

Accurate laser altimeter elevation profiles contribute to a 
number of Solid Earth science and applied objectives. A 
primary contribution is the independent characterization of 
systematic and random elevation errors in Digital Elevation 
Models (DEMs) produced by photogrammetric and 
Interferometric SAR techniques. This kind of quality 
assessment enables DEMs to be used quantitatively for 
purposes for which their accuracy is appropriate. In addition 
to evaluating DEM accuracy, laser altimeter profiles can be 
used in the correction of systematic errors in DEMs, im- 
proving their utility for detection of elevation change 
observed by differencing DEMs obtained at different times. 
NASA’s Deformation, Ecosystem Structure and Dynamics 
of Ice (DESDynl) mission, scheduled for launch later this 
decade, will provide globally distributed laser altimeter 
profiles well suited for these purposes. It will 
simultaneously acquire multiple profiles with 25 m 
diameter, nearly-contiguous laser footprints along the 
profiles. Although NASA’s Ice, Cloud and land Elevation 
Satellite (ICESat) has less dense sampling (a single profile 
with ~50 m footprints spaced 175 m apart), it illustrates the 
capabilities for accuracy assessment and control of DEMs 
that will be significantly expanded by DESDynl. 

We have previously used ICESat data to evaluate the 
accuracy of the 90 m resolution, near-global DEM produced 
by the Shuttle Radar Topography Mission (SRTM) [4, 5, 6] 
and the more recent global ASTER DEM [7]. The 
Geoscience Laser Altimeter System (GLAS) aboard ICESat 
acquired data from 2003 to 2009, during month-long 
observant periods that were conducted three, and later two, 
times a year. Figure 1 shows an example assessment of 
SRTM accuracy using an ICESat profile. The negative- 
elevation differences reveal that the SRTM DEM is biased 
high relative to an absolute datum by several meters, on 
average, across western Australia. In addition, the along- 
profile variations reveal undulating elevation errors in the 
SRTM DEM at the 100s of kilometer length scale and ~5 m 
amplitude. ICESat data documents the spatial structure of 
these undulations, which can then be removed from the 
DEM. The comprehensive, near global (± 86°) ICESat 


coverage across continents enables correction of long- 
wavelength DEM errors not previously made possible by 
other means. Here we present results for along-profile and 
gridded elevation differences between ICESat geodetic 
control data and the SRTM DEM for Australia, 
documenting the spatial pattern and magnitude of elevation 
errors in that DEM. We do this using each ICESat 
observation period separately, documenting that the ICESat 
results are highly reproducible. 



Figure 1. Elevation measurements (top right) along an ICESat 
profile in western Australia (left, red line on MODIS image) and 
differences with the SRTM DEM (bottom right) for individual 
ICESat laser footprints. 

2. DEVELOPMENT OF A GLOBAL ICESAT 
GEODETIC CONTROL DATABASE 

Using the ICESat Land/Canopy Elevation data product 
(GLAM), we are producing a global set of GCPs in a 
project supported by NASA’s Earth Surface and Interior 
Program. We apply stringent editing criteria in order to 
produce the highest quality ground control. For this paper, 
we have used the latest public release of the ICESat 91 -day 
repeat track data (Release 31). Like the DESDynl lidar will 
do, the GLAS instrument records a waveform, the height 
distribution of energy reflected from illuminated vegetation 
and ground surfaces in the laser footprint. In addition to 
providing the latitude, longitude and derived elevations for 
each ICESat waveform referenced to the Topex/Posseidon 
(T/P) ellipsoid, the GLAM product also includes the SRTM 
DEM elevation at the footprint location. The publically 
distributed finished SRTM DEM contains orthometric 
elevations produced using the EGM96 geoid [8]. On the 
GLAM Release 31 product that we are using the SRTM 
elevations have been converted to ellipsoidal values, but by 
incorrectly applying the more recent EGM08 geoid. We 
therefore correct the SRTM ellipsoidal heights in the 
products using the SRTM reference geoid, EGM96. We then 
compute the elevation differences between the SRTM radar 
phase center elevation and four elevations derived from the 
ICESat waveform: the highest observed surface (waveform 
start), the waveform centroid (average elevation), the lowest 
observed surface (waveform end) and, where possible, the 
ground elevation inferred from the lowest Gaussian fit 
derived from waveform modeling. We also obtain a land 
cover classification for each laser footprint from the MER1S 


Globcover land cover regional product (5 1 possible classes), 
derived by an automatic and regionally-tuned classification 
of a MERIS FR time series [9]. The MERIS product is 
generated at 300 m resolution using data from the period 
December 2004 - June 2006. 

We apply stringent editing criteria to yield a high 
quality GCP database. We exclude ICESat data identified as 
returns from water (MERIS land cover = 210). We also 
exclude ICESat returns inferred to be from clouds using a 
threshold 50 m above the SRTM surface at the laser 
footprint location. The cloud flag in the ICESat products 
indicating a clear atmosphere proved too severe for editing 
and was not reliable for periods when the transmit laser 
energy was low. One of the main criteria for obtaining the 
most accurate ground control data is selection of waveforms 
from relatively non-vegetated areas with low relief. We do 
so by limiting the width of the selected waveforms (signal 
start to signal end) to 5 m. For these narrow waveforms the 
elevations for the waveform centroid and the ground lowest 
Gaussian fit essentially coincide. 

We also apply editing based on instrumental 
parameters. Weak returns are excluded by only using returns 
with maximum amplitudes greater than 0.15 Volts. Where 
the return energy exceeds the receiver’s dynamic range, the 
GLAS waveforms become saturated. Saturation broadens, 
distorts and truncates the received waveforms. Although 
laboratory calibrations developed as a function of receiver 
gain and observed received energy exist to correct 
elevations and receive energies, we only keep non-saturated 
to minimally saturated data where the saturation index (the 
number of waveform bins with an amplitude greater than the 
saturation threshold) is two or less. The accuracy of the 
ICESat data is degraded with increasing incidence angle 
between the laser beam vector and the normal to the surface 
slope, causing waveform broadening. To minimize this error 
source, we exclude data acquired when the laser beam was 
pointed off from nadir by more than 1°. Rigorous analysis 
has shown that for low relief locations the ICESat data meet 
the accuracy requirements of 6 m horizontal and 10 cm 
vertical [S. Luthcke, pers. comm.]. Because we use 
stringent editing criteria we expect that our GCPs are of 
equivalent accuracy. 

Table 1 presents elevation difference between our 
ICESat GCPs and SRTM, computed as ICESat minus 
SRTM, for Australia. Results are given for ICESat 
observation periods acquired with the 2 n ° (L2) and 3 rd (L3) 
laser transmitters, with their start and end dates and average 
laser transmit energies indicated (for each laser the transmit 
energy started high, causing substantial saturation, and then 
dropped significantly during the course of the mission). 

Figure 2 illustrates ICESat minus SRTM difference 
histograms for ICESat’s highest, centroid and lowest 
elevations using a representative observation period (L3E), 
showing well-defined normal distributions. Like the other 



periods, this period shows an SRTM elevation bias with 
respect to the ICESat centroid of ~ -2 m. 
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Tabic 1. ICESat- SRTM Elevation Difference Statistics 

D = Differences (H: Highest; C: Centroid; L: Lowest; G: ground) 

N* = Number of returns classified as be from the Earths’s surface 
and the percent retained after editing 
N = Number of returns after editing 

M =Mean elevation difference: S = standard deviation; R = RMSE 
E (fl)** = Average transmit energy corrected for receiver field of view 
(FOV) shadowing 



Figure 2. Elevation difference histograms for ICESat highest, 
centroid and lowest minus SRTM for the L3E laser observation 
period (See Table 1 for complete statistics). 

Examination of the centroid differences for all laser 
periods (Figure 3) shows very consistent mean elevation 
difference results, a demonstration of ICESat’s highly 
accurate and reproducible absolute elevations. There is a 
slightly decreasing trend with laser energy decay, especially 
for Laser 2 (L2), potentially related to the large number of 
edited saturated returns during early, high energy periods. 
This could cause a bias by preferentially removing areas of 
high surface reflectivity. The origin of this L2 drift and the 
associated increase in standard deviation requires further 
investigation. 
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Figure 3. Mean ICESat centroid minus SRTM elevation 
differences and standard deviations for all ICESat observation 
periods (See Table 1 for complete statistics). 

3. DERIVED CORRECTION SURFACES FOR 
AUSTRALIA 

After selecting the ground control database for each 
ICESat period, we obtain a representation of the undulations 
in ICESat centroid minus SRTM differences in two different 
ways. First, we apply a sliding boxcar filter along every 
edited ICESat profile at 1° intervals (Figure 4, top). 
Alternatively, we divide the region in 1° by 1° tiles, and 
compute the mean differences and their standard deviations 
for the points within each tile (Figure 4, lower 3 panels). 
Both the along-profile and gridded representations highlight 
long wavelength elevation-error undulations in the SRTM 
DEM of several meters magnitude. 


4. CONCLUDING REMARKS 
AND APPLICATION TO THE DSDYNI MISSION 
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Figure 4. Mean ICESat centroid - SRTM differences for the L3E 
observation period. Differences computed using a 1° along track 
sliding box-car Filter are shown on the top, and 1° gridded averages 
and their standard deviations for the indicated number of returns in 
a cell are shown in the three lower panels. Colors indicate ranges 
shown by the color bars. 


Profiling spaceflight laser altimeters make important 
contributions to the solid Earth Science by providing 
accurate, global, consistently referenced, geodetic elevation 
data. Expanding upon the Australia work described here, 
our ICESat GCP dataset is contributing to a global, 
coordinated and integrated DEM database produced from 
different sources, embedded into a consistent, high 
accuracy, and long term stable geodetic reference frame. It 
will be particularly useful in northern and southern latitudes 
above and below ±60°, where high-resolution topographic 
data assets are not available and topographic control is 
scarce. Methodologies developed to use ICESat data for 
global geodetic control purposes are a pathfinder for similar 
use of data to be produced by the lidar component of the 
DESDynl mission. With substantially improved sampling 
as compared to ICESat, DESDynl will provide a more 
comprehensive set of global GCPs from it’s multi-beam, 
higher-resolution elevation profiles. 
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